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Используя пространственную геомеханическую модель деформирования и разрушения горных пород в зоне влияния горной 
выработки, разработанную на основе метода конечных элементов, проведено исследование влияния движущегося очистного 
забоя на характер изменения коэффициента концентрации вертикальных напряжений впереди очистного забоя, зависания и 
циклического обрушения подработанных пород кровли угольных пластов. 


Введение 

В последнее десятилетие резко возросла интен- 
сивность подземной разработки угольных место- 
рождений. За этот период нагрузка на очистной 
комплексно-механизированный забой возросла в 
3...4 раза, а длина лавы увеличилась от 120 до 300 м. 
Закономерности геомеханических процессов для 
таких забоев изучены недостаточно, что подтверж- 
дается большой аварийностью забоев и высоким 
уровнем группового травматизма (шахты "Зыряно- 
вская", "Тайжина" в Кузбассе и др.). 

Изучением геомеханических процессов на 
угольных шахтах в окрестности очистных забоев за- 
нимались многие ученые: Ф.А. Белаенко, А.А. Бо- 
рисов, Г.Н. Кузнецов, В.Д. Слесарев, П.М. Цимба- 
ревич, И.Л.Черняк и др. Общие закономерности 
сдвижения подрабатываемого массива горных по- 
род неоднократно описаны в литературе [ 1 -5 и др. ] . 

Механизм разрушений и обрушений толщи по- 
род кровли разработан Г.Н. Кузнецовым [1]. Он до- 
казал, что обрушение слоев непосредственной 
кровли в выработанном пространстве будет иметь 
место при условии, если толщина слоев кровли 
меньше мощности пласта. Характер первого обру- 


шения основной кровли в срединном сечении вы- 
работанного пространства, ее взаимодействие с не- 
посредственной кровлей и характер опорного дав- 
ления на краевые зоны разрабатываемого пласта 
схематично показаны на рис. 1. 

Как отмечается в [2] породы над зоной интен- 
сивного разрыхления деформируются и разруша- 
ются подобно тонким плитам, в режиме первого 
обрушения - подобно плите, защемленной на пря- 
моугольном опорном контуре, в режиме устано- 
вившегося обрушения - подобно консольной пли- 
те, защемленной по трем сторонам опорного кон- 
тура. 

Общий характер распределения опорного дав- 
ления в плоскости пологого пласта, рассчитанный 
по результатам физического моделирования на 
объемных моделях, проверенный и многократно 
подтвержденный шахтными и лабораторными ис- 
следованиями описан А. А. Борисовым [3, 4]. По 
характеру изменений опорного давления он выде- 
ляет следующие зоны: динамических проявлений 
опорного давления, затухания динамических про- 
явлений опорного давления, статического состоя- 
ния опорного давления. 


39 


Известия Томского политехнического университета. 2005. Т. 308. № 1 



Рис. 1 . Характер первого обрушения основной кровли в среднем сечении ( по простиранию) и распределение опорного давле- 
ния на краевые зоны разрабатываемого пласта по Г.Н. Кузнецову [1] 


Анализ закономерностей проявления горного 
давления аналитическими методами выполнен 
Г.Н. Кузнецовым, А.М. Линьковым, И.М. Петухо- 
вым, К.В. Руппенейтом, Г.Н. Савиным и другими 
авторами. В аналитических методах исследования 
вопросов механики горных работ для расчета нап- 
ряженного состояния массива в зонах влияния 
подземных выработок используются положения 
теории упругости, пластичности, ползучести. Учет 
влияния неоднородности состава пород междуп- 
ластий на напряженное состояние толщи и уголь- 
ных пластов выполнялся по методу эквивалентно- 
го слоя, широко применяемого в механике грунтов. 

Разработанные в настоящее время аналитичес- 
кие методы расчета напряжений в зонах влияния 
подземных выработок имеют ограниченную об- 
ласть применения и используются, главным обра- 
зом, на стадии прогнозной оценки горно-геологи- 
ческих и горнотехнических условий отработки, для 
выбора способов и средств управления горным 
давлением. 

Авторами [5] были проведены исследования ха- 
рактерных зон и параметров сдвижения массива 
под действием очистных работ на моделях из экви- 
валентных материалов. Измерения тензодатчика- 
ми вертикальных напряжений в подработанном 
массиве показали, что в 30. ..40 м впереди забоя 
напряжения возрастают в 1,5. ..2,0 раза, затем над 
забоем, или с некоторым опережением забоя, про- 
исходит снижение уровня напряжений практичес- 
ки до нуля. В дальнейшем происходит восстанов- 
ление напряжений. Полностью до исходного уров- 
ня напряжения восстанавливаются на расстоянии 
50. ..60 м за забоем. 

Следует отметить, что одним из недостатков су- 
ществующих методов прогнозирования напряжен- 
но-деформированного состояния углепородного 
массива при проведении горных выработок являет- 
ся решение двумерной задачи, что не обеспечивает 
корректного результата на сопряжении очистного 
забоя и угольного массива. 

Кроме того, сложность расчетов напряженного 
состояния толщи пород при разработке свит плас- 
тов обусловлена необходимостью учета взаимного 


влияния выработок на различных пластах и неод- 
нородностью состава углепородного массива. В ре- 
альных условиях учитывать влияние всех слоев гор- 
ных пород на напряженное состояние толщи весь- 
ма трудоемко и сложно из-за недостаточной на- 
дежности определения механических характерис- 
тик на контактах слоев. 

В этой связи научные исследования геомехани- 
ческих процессов в окрестности высокопроизводи- 
тельного комплексно-механизированного забоя 
путем решения пространственной задачи геомеха- 
ники с использованием численных методов явля- 
ются актуальными. 

Объект и методы исследований 

Для сравнения фактических данных с числен- 
ными результатами количественного прогнозиро- 
вания геомеханических процессов в углепородном 
массиве при движении очистного забоя в качестве 
базового варианта объекта исследований выбран 
выемочный участок 26-26 шахты "Абашевская" в 
Кузбассе. 

При расчете параметров напряженно-деформи- 
рованного состояния углепородного массива в зоне 
влияния горной выработки использовалась простра- 
нственная геомеханическая модель деформирования 
и разрушения горных пород, разработанная автора- 
ми на основе метода конечных элементов [6]. Прог- 
нозирование геомеханических процессов в углепо- 
родном массиве выполнялось на основе анализа 
пространственного напряженно-деформированного 
состояния и использования обобщенного критерия 
прочности Кулона-Мизеса. При моделировании ме- 
ханизма обрушения горных пород осуществлялось 
последовательное накопление повреждений с уче- 
том функции времени. Для имитации движения 
очистного забоя выполнялся циклический сдвиг 
системы координат по мере продвижения забоя. 

Учитывая сложность поставленной задачи и ог- 
раничения, связанные с техническими возможнос- 
тями вычислительной техники, проводились ис- 
следования характера изменения коэффициента 
концентрации вертикальных напряжений впереди 
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Рис. 2. Схема движения очистного забоя, ограниченного по восстанию и падению угольным массивом на: а) 10; 6) 40 м 


очистного забоя, зависания и циклического обру- 
шения подработанных пород кровли при движении 
длинного очистного забоя, ограниченного по вос- 
станию и падению угольным массивом. Имитация 
движения очистного забоя от монтажной камеры 
выполнялась с шагом 10 м (рис. 2, а, б). 

На каждом шаге моделирования вычислялись 
горизонтальные и вертикальные смещения, дефор- 
мации, напряжения, определялись коэффициент 
концентрации вертикальных напряжений и оста- 
точная прочность угля и пород в зоне влияния вы- 
работанного пространства. Наиболее представи- 
тельными оказались графики изменения коэффи- 
циента концентрации вертикальных напряжений в 
угольном массиве впереди очистного забоя, а так- 
же формы и размеры зон зависания и обрушения 
подработанных пород кровли. 

Исследование влияния движущегося забоя 
на характер зависания и циклического обрушения 
подработанных пород кровли угольных пластов 

В результате проведенных исследований уста- 
новлено, что при движении очистного забоя зона 
обрушения пород кровли имеет ассиметричную 
форму. По мере движения забоя от монтажной ка- 
меры происходит последовательное обрушение по- 
род непосредственной кровли в середине вырабо- 
танного пространства. Высота зоны обрушения 
увеличивается пропорционально расстоянию от 
монтажной камеры до забоя. Наиболее интенсивно 
разрушение породы происходит вблизи движуще- 
гося очистного забоя (рис. 3). 



Рис. 3. Схема формирования зоны обрушенных пород под- 
работанной кровли 

Над угловыми участками выработанного прост- 
ранства происходит зависание пород непосред- 
ственной кровли, что приводит к увеличению нап- 
ряжений в угольном массиве. 

При наличии в подрабатываемой толще уголь- 
ных пластов разрушение угля в них опережает раз- 
рушение пород непосредственной кровли. Над вы- 
работанным пространством формируется порода- 
мост, мощностью равной расстоянию от кровли от- 
рабатываемого пласта до подрабатываемого пласта, 
т.е. происходит расслоение подрабатываемой тол- 
щи и формирование двух породных плит (рис. 4). 

Наличие таких пород-мостов над выработан- 
ным пространством приводит к росту напряжений 
в отрабатываемом пласте, увеличению длины кон- 
соли при движении забоя и динамическому обру- 
шению породы-моста при достижении предельно- 
го значения длины консоли. 
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Рис. 4. Схема формирования зависающей консоли подра- 
ботанных пород кровли 


Динамическое обрушение пород кровли сопро- 
вождается заполнением выработанного простран- 
ства, вытеснением метано-воздушной смеси и 
угольной пыли в горные выработки и, при наличии 
источника открытого огня, подземным взрывом. 


Исследование влияния движущегося забоя 
на характер изменения коэффициента 
концентрации вертикальных напряжений 
впереди очистного забоя 

По мере движения от монтажной камеры проис- 
ходит почти линейное увеличение коэффициента 
концентрации вертикальных напряжений (рис. 5). 
Прирост вертикальных напряжений происходит за 
счет увеличения веса зависающих пород кровли. 

При движении очистного забоя одновременно с 
зависанием пород основной кровли и налегающей 
толщи происходит обрушение пород непосред- 
ственной кровли, что приводит к частичному сни- 
жению веса зависающих подработанных пород. 
Это обрушение пород непосредственной кровли 
зафиксировано на графиках на рис. 5 в виде сниже- 
ния коэффициента концентрации вертикальных 
напряжений после каждого шага подвигания забоя. 

Как следует из графиков на рис. 5, снижение ко- 
эффициента концентрации вертикальных напря- 
жений в среднем сечении лавы (точка В на рис. 2, а) 
приводит к его возрастанию над угловыми участка- 
ми лавы (точки А и С на рис. 2, а). 

В процессе движения очистного забоя на рас- 
стоянии равном 0,7 длины лавы резко возросла 


■■ — сопряжение со штреком (точка А) середина лавы (точка В) 
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Рис. 5. Г оафики изменения коэффициента концентрации вертикальных напряжений при движении очистного забоя 
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нагрузка на угольный пласт в средней части лавы. 
На следующем шаге произошло первичное обру- 
шение пород основной кровли, что привело к разг- 
рузке и снижению коэффициента концентрации 
вертикальных напряжений до 0,87 в призабойной 
части угольного массива в середине лавы (рис. 5). 

После частичного обрушения пород основной 
кровли в середине выработанного пространства 
произошло перераспределение веса необрушенных 
пород на краевые участки угольного пласта, т.е. за 
счет взаимного влияния линии забоя и границ вы- 
емочного участка в угловых частях выработанного 
пространства кровля не обрушилась, и резко воз- 
росло горное давление. Это привело к увеличению 
коэффициента концентрации вертикальных нап- 
ряжений в пласте до 3,22 вблизи сопряжения вен- 
тиляционного штрека и линии забоя (рис. 5). 

При следующем подвигании забоя на 10 м про- 
изошло полное обрушение пород кровли в преде- 
лах всего отработанного контура. Коэффициент 
концентрации вертикальных напряжений в уголь- 
ном пласте вблизи сопряжений с горными выра- 
ботками составил 1,48 (рис. 5). 

В процессе дальнейшего подвигания лавы сно- 
ва начала формироваться консоль зависания пород 
кровли. Коэффициент концентрации вертикаль- 
ных напряжений в угольном массиве увеличивает- 
ся как в середине лавы, так и на верхнем и нижнем 
сопряжениях лавы с горными выработками. При 
движении забоя происходит частичное обрушение 
пород непосредственной кровли, что приводит к 
снижению коэффициента концентрации верти- 
кальных напряжений. 

Через 30 м после первичного обрушения коэф- 
фициент концентрации вертикальных напряжений 
увеличился примерно в 1,3 раза, и произошло вто- 
ричное обрушение. Вторичному обрушению соот- 
ветствует резкое уменьшение коэффициента кон- 
центрации вертикальных напряжений в 1,4 раза 
(рис. 5). Длина консоли вторичного зависания сос- 
тавила 0,33 от длины первичного шага обрушения. 


Ввиду высокой аварийности горно-шахтного 
оборудования возможны остановки очистных за- 
боев продолжительностью до нескольких месяцев. 
В этом случае вследствие высоких напряжений в 
угольном пласте вблизи от забоя происходит их ре- 
лаксация, т.е. снижение коэффициента концентра- 
ции вертикальных напряжений (рис. 6). 

Релаксация напряжений сопровождается 
уменьшением прочности угля, что приводит к рос- 
ту смещений пород кровли над очистным забоем. 
При этом вследствие деформации пород кровли и 
пучения пород почвы происходит конвергенция 
перекрытия и основания секций крепи, "посадка" 
стоек механизированной крепи насухо и возникно- 
вение аварийной ситуации. 

По результатам расчетов также установлено, что 
при остановке очистного забоя происходит увели- 
чение высоты зоны обрушения пород кровли, в том 
числе с образованием трещин и заколов впереди 
очистного забоя. 

Таким образом, полученные результаты удов- 
летворительно описывают циклический характер 
зависания и обрушения подработанных пород 
кровли при движении очистного забоя. Отношение 
вторичного и первичного шагов обрушения попа- 
дают в диапазон 0,25. ..0,5, установленный экспери- 
ментально [1-5]. 

Изменение вычисленного коэффициента кон- 
центрации вертикальных напряжений в диапазоне 

1.3. . .2. 15 соответствует экспериментальным дан- 
ным, диапазон значений которых составляет 

1.4.. .2. 29 по замерам напряжений на пластах впере- 
ди очистных и подготовительных выработок гид- 
равлическими датчиками, выполняемым на протя- 
жении ряда лет в ВостНИИ [7]. 

Разработанная авторами пространственная гео- 
механическая модель деформирования и разруше- 
ния горных пород позволяет прогнозировать пре- 
дельное время остановки очистного забоя в случае 
аварийной ситуации с учетом раздвижности сек- 
ций механизированной крепи. 
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Рис. 6. Г оафики изменения коэффициента концентрации вертикальных напряжений при остановке очистного забоя 


43 





Известия Томского политехнического университета. 2005. Т. 308. № 1 


СПИСОК ЛИТЕРАТУРЫ 

1. Кузнецов Г.Н. Механические свойства горных пород. — М.: Уг- 
летехиздат, 1947. — 180 с. 

2. Черняк И.Л., Ярунин С.А. Управление состоянием массива 
горных пород. — М.: Недра, 1995. —395 с. 

3. Борисов А.А. Механика горных пород и массивов. — М.: Нед- 
ра, 1980. —360 с. 

4. Борисов А.А. Расчеты горного давления в лавах пологих плас- 
тов. — М.: Недра, 1964. —278 с. 

5. Управление горным давлением при разработке пологих плас- 
тов с труднообрушаемой кровлей на шахтах Кузбасса / 


С. И. Калинин, А.Ф. Лютенко, В.П. Егоров, С.Г. Дьяконов. — 
Кемерово: Кемеровское книжное изд-во, 1991. —247 с. 

6. Павлова Л.Д., Фрянов В.Н. Моделирование циклического ха- 
рактера обрушения горных пород при проведении выработки с 
последовательным накоплением повреждений // Известия 
Томского политехнического университета. — 2004. — Т. 307. — 
№2. -С. 76-79. 

7. Чернов О.И., Мурашев В. И., Шлиомовичус Я.Г. Распределе- 
ние напряжений вокруг выработок с учетом псевдопластичес- 
ких свойств угольного массива // Проблемы механики горных 
пород: Матер. Всес. научной конф. по механике горных пород. 
— Новосибирск: Наука, 1971. — С. 684-689. 


ѵліе™ ю-чі 4/™ пкі 




